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Abstract 

A description of neutral and multiply charged fullerenes is proposed based on a stabilized jel- 
lium (structureless pseudopotential) approximation for the ionic background and the local density 
approximation for the a and vr valence electrons. A recently developed shell-correction method is 
used to calculate total energies and properties of both the neutral and multiply charged anionic 
and cationic fullerenes. The effect of the icosahedral symmetry is included perturbatively. The 
calculated single-particle energy level spectrum of Cgo is in good correspondence with experimen- 
tally measured ones and previous self-consistent local-density-approximation calculations. For the 
multiply charged fullerenes, we calculate microscopically the charging energies of Cg^ for up to 
X = 12 excess charges. A semiclassical interpretation of these results is developed, which views the 
fullerenes as Coulomb islands possessing a classical capacitance. The calculated values for the first 
ionization potential and the first electron affinity agree well with the experimental ones. For the 
second and third ionization potentials, there exist substantial discrepancies in the experimental 
measurements. Our calculations support the results from charge transfer bracketing experiments 
and from direct ionization experiments through electron impact. The doubly charged negative ion 
is found to be a very long-lived metastable species, in agreement with observations. 

PACS numbers: 
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INTRODUCTION 



Charging of macroscopic metal spheres is an old subject with scientific accounts dating 
back to Coulomb, Faraday, and others Recently, issues related to electrical charging 
emerged in connection with quantal nanostructures which are of interest to diverse areas 
of condensed-matter and atomic, molecular and cluster physics, such as nanofabricated 
semiconductor devices known as Coulomb islands 3, sl, metal microclusters {4, Is], O, 3], 



and fullerenes 
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15, 
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13, 



electron attachment to molecular [8| , [9| and alkali-halide [l^ clusters, and carbon clusters 



18|- 



We report on a unified theoretical study of the energetics, stability, ionization potentials 
and electron affinities of neutral and charged carbon fullerenes, Cg^, using local density 
functional theory (LDA) 19| with a stabilized-iellmm background \2Q]. The present approach 
is an adaptation to the case of fullerenes of a shell-correction method (SCM), developed by 
us previously for metal microclusters . The main elements of this adaptation consist of 
employment of the stabilized jellium instead of the usual jellium background, and the use 
of a generalized electron-density profile (appropriate for description of the fullerene cage). 
Furthermore, the point-group icosahedral symmetry of the Cgo cage is introduced via a 
perturbative treatment. 

The electronic structure of fullerenes has been investigated extensively using various 
methods ranging from ab initio all-electron quantum chemical Hartree-Fock calculations, and 



self-consistent Kohn-Sham LDA calculations in conjunction with non 



oca. 



pseudopotentials 
21I . 122( 1 and particle-on- 

231]. Our generalized shell-correction method (GSCM) which combines 



or a jellium background for the ions, to simplified free-electron 
a-sphere models 



elements of the LDA methodology with simplifications circumventing self-consistent solution 
of the Kohn-Sham (KS) equations represents a substantial improvement over certain previous 
simplified methods, yielding results in quantitative agreement with self-consistent KS-LDA 



24 



25 



26|). and experimental data. 



ab initio pseudopotential calculations (see, e.g.. 

Additionally, our method permits us to investigate, within the local density approxima- 
tion, the important class of multiply anionic fullerene systems, where the familiar KS-LDA 
is known to fail 0, Until now, in all cases but one 27], properties of anionic fullerenes 
(whose energetics is relevant not only for understanding the properties of free molecules, but 
also those of fullerite intercalation compounds 29] and complexes js^) have been treated 
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within the Hartree-Fock approximation 



cal electrostatic model 
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31| or estimated with a simplified classi- 



29|, associated with the particle-on-a-sphere model. Both these 



methods, however, omit essential theoretical ingredients, namely the correlation energy in 
the case of the former, and the shell structure and appropriate radial and angular density 
distribution of the electronic charge in the case of the latter. The recent study by Pederson 
and Quong {2^, using an LDA all-electron full-potential Gaussian-orbital basis, was per- 
formed within the ansatz of a localized basis expansion of the Kohn-Sham orbitals, yielding 
results in good agreement with those obtained by us using the GSCM. 



LDA THEORETICAL METHOD 



Stabilized jellium approximation 



FuUerenes and related carbon structures have been extensively investigated using ab initio 
local-density- functional methods and self-consistent solutions of the Kohn-Sham (KS) equa- 



tions 
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321 ]. For metal clusters, replacing the ionic cores with a uniform jellium background 



was found to describe well their properties within the KS-LDA method (see references in 
Refs. 0, ^). Motivated by these results, several attempts to apply the jellium naodel in 



conjunction with LDA to investigations of fuUerenes have appeared recently 
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35|. 



Our approach differs from the earlier ones in several aspects and, in particular, in the adap- 
tation to the case of finite systems of the stabilized-jellium (or structureless pseudopotential) 



energy density functional (see eq. ([T]) below and Ref. |20l|).|45| 

An important shortcoming of the standard jellium approximation for fuUerenes (and other 
systems with high density , i.e., small r^) results from a well-known property of the jellium 
at high electronic densities, namely that the jellium is unstable and yields negative surface- 
energy contribution to the total energy 2^ , as well as unreliable values for the total energy. 
These inadequacies of the standard jellium model can be rectified by pseudopotential cor- 
rections. A modified-jellium approach which incorporates such pseudopotential corrections 
and is particularly suited for our purposes here, is the structureless pseudopotential model 
or stabilized jellium approximation developed in Ref. 2o| . 

In the stabilized jellium, the total energy Epseudo, as a functional of the electron density 



3 



p(r), is given by the expression 

Epseudo[p, P+] = Ejeii[p,p+] + {Sv)ws J p{r)U{r)dT -e j p+{r)dr , (1) 

where by definition the function U{r) equals unity inside, but vanishes, outside the jelhum 
volume. P4. is the density of the positive jellium background (which for the case of Ceo 
is taken as a spherical shell, of a certain width 2d, centered at 6.7 a.u. ). Epseudo in eq. 
([1]) is the standard jellium- model total energy, Ejeii, modified by two corrections. The 
first correction adds the effect of an average (i.e., averaged over the volume of a Wigner- 
Seitz cell) difference potential, {6v)ws^{^), which acts on the electrons in addition to the 
standard jellium attraction and is due to the atomic pseudopotentials (in this work, we 



use the Ashcroft empty-core pseudopotential, specified by a core radius r^, as in Ref. |20l|). 
The second correction subtracts from the jellium energy functional the spurious electrostatic 
self-repulsion of the positive background within each cell; this term makes no contribution 
to the effective electronic potential. 



Following Ref. [20|], the bulk stability condition (eq. (25) in Ref. [20|]) determines the 
value of the pseudopotential core radius function of the bulk Wigner-Seitz radius r,. 

Consequently, the difference potential can be expressed solely as a function of as follows 



(energies in Ry, distances in a.u. 



where Ec is the per particle electron-gas correlation energy (in our calculation, we use the 



3))- 



Gunnarsson-Lundqvist exchange and correlation energy functionals (see Refs. 

The electrostatic self-energy, e, per unit charge of the uniform positive jellium is given by 

e = 6v^/y5rs, (3) 

where v is the valence of the atoms {v = 4 for carbon). 

Shell-correction method 

Besides the familiar KS-LDA approach, an alternative LDA method, which has been used 



in studies of metal clusters 



36|, is based on an extended Thomas Fermi (ETF) vari- 



ational procedure using a parametrized density profile p(r;{7j}), with {•ji} as variational 



parameters 5|]. In the case of metal clusters, the energy density functional that is varia- 
tionally minimized consists of a kinetic-energy functional T[p], including terms up to 4th 
order in the density gradients (see Ref. 5(]), and of potential terms according to the standard 
jellium-LDA functional. The effective single-particle potentials and associated single-particle 
energy spectra obtained by this method for metal clusters provide a good approximation 
to the corresponding ones obtained from KS-LDA calculations 0, ^, 



and have been exten- 



sively used 



36l | in studies of the optical properties of metal clusters. To apply the ETF-LDA 
method to carbon fuUerenes, we generalize it by employing potential terms according to the 
stabihzed-jellium functional ([I]). 

Another required generalization consists in employing a parametrized electron-density 
profile that accounts for the hollow cage-structure of the fuUerenes. Such a density profile 
is provided by the following adaptation of a generalization of an inverse Thomas-Fermi 



distribution, used earlier in the context of nuclear physics 37|, i.e. 



° \cosh[wi^o/ai,o] + cosh[(r - R)/ai^o\J 
where R = 6.7 a.u. is the radius of the fullerene cage, w, a, and 7 are variables to be 
determined by the ETF-LDA minimization. For R = and large values of w/a, expression 
(jlj) approaches the more familiar inverse Thomas-Fermi distribution, with w the width, a 
the diffuseness and 7 the asymmetry of the profile around r = w. There are a total of six 
parameters to be determined, since the indices {i, o) stand for the regions inside (r < R) and 
outside (r > R) the fullerene cage. Fi^o = (cosh[wj o/aj^o] + 1)/ sinh[wj^o/Q;j J is a constant 
guaranteeing that the two parts of the curve join smoothly at r = i?. The density @ peaks 
at r = and then falls towards smaller values both inside and outside the cage (see top 
panel of Fig. 1). 

However, apart from the effective potentials, other properties of metal clusters - and by 
analogy of fuUerenes - determined via the ETF-LDA method, such as ionization potentials 
(IPs), electron affinities (EAs), and total energies, compare with the self-consistently calcu- 
lated KS-LDA values only in an average sense, i.e., they do not exhibit shell effects. In our 

E resent study, shell effects are incorporated in the GSCM according to the procedure of Refs. 
in n 
y, 15| . As elaborated in Ref. [5| , shell effects in the total energy are contained, to first order 

in (5pKS {^Pks = Pks — P, where p is the ETF-LDA optimized density specified above and pks 

is the self-consistent KS density), in the sum ^^^i, where Ei are the single-particle energies 
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FIG. 1: Bottom panel: The stabilized-jellium LDA potential obtained by the ETF method 
for the neutral Ceo molecule. The Wigner-Seitz radius for the jellium bacground is 1.23 a.u. 
Note the asymmetry of the potential about the minimum. The associated difference potential 
{Sv)ws = -9.61 eV. 

Top panel: Solid line: Radial density of the positive jellium background. Dashed line: ETF 
electronic density. Note its asymmetry about the maximum. Thick solid line: The difference (mul- 
tiplied by 10) of electronic ETF densities between CgQ and Cgo- It illustrates that the excess charge 
accumulates in the outer perimeter of the total electronic density. All densities are normalized to 
the density of the positive jellium background. 

of the ETF effective potential, and the sum extends over the occupied levels. Consequently, 
we replace the kinetic-energy T in the ETF-LDA functional by 



where V{r] p{r)) is the effective potential produced by the ETF method [4, |5| (as an example, 
see lower panel Fig. 1 for the ETF potential associated with the neutral Ceo)- As a result, 
the total energy. Est, including the shell correction, Ai^sh = Tg^^ — T, is given by -Esh[p] = 




(5) 
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— T + -E'[p], where E is the ETF-LDA energy-density functional. 

After some rearrangenments, the shell-corrected total energy -Esh[p] in the GSCM can be 
written in functional form as follows 

Est[p\ = J2^^- j + Kc(r)| p(r)dr + j S.,[p{v)]dv + Ei-e j p+{v)dv , (6) 

where Vh and V^c are the Hartree and exchange-correlation electronic potentials, £^c is the 
exchange- correlation energy density functional, and Ej is the energy of the positive jellium 
background. The specific way of writing the functional E^h above was chosen so that its 
similarity in form to the Harris functional 38|] is evident. We note that our method differs 



from that approach in that the optimization of the input density is achieved by us through 
a variational ETF method, which does not require a step-by-step matrix diagonalization. 
While our focus in the previous papers 0, was on jellium models for metal clusters, the 
very good agreement between our results and those obtained via self-consistent KS-LDA 
jellium calculations 0, suggested that it would be worthwhile to explore the application 
of our SCM to more general electronic systems extending beyond metal clusters. The present 
study is an example of such an application of the SCM to fuUerenes. 

Icosahedral splitting 

Heretofore, the point-group icosahedral symmetry of Cgo was not considered, since the 
molecule was treated as a spherically symmetric cage. This is a reasonable zeroth-order 



approximation as noticed by several authors 22|, |2J, |35|, l39|] . However, considerable improve- 



ment is achieved when the effects of the point-group icosahedral symmetry are considered 
as a next-order correction (mainly the lifting of the angular momentum degeneracies). 

The method of introducing the icosahedral splittings is that of the crystal field theory. 
Thus, we will use the fact that the bare electrostatic potential from the ionic cores, considered 



as point charges, acting upon an electron, obeys the well-known expansion theorem 401] 

U{r) = -ve' J2 u— 77 = - E E '^iir)CrYr{0, 4>) , (7) 

i 1=0 m=—l 

where the angular coefficients are given through the angular coordinates 6'j, (pi of the 
carbon atomic cores, namely, 

cr = Y.^r{o.A^) , (8) 
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FIG. 2: (a) The single-particle levels of the ETF-LDA potential for Cgo shown in Fig. 1. Because 
of the spherical symmetry, they are characterized by the two principle quantum numbers n,. and I, 
where is the number of radial nodes and I the angular momentum. They are grouped in three 
bands labeled a {rir = 0), vr {rir = 1), and 6 [rir = 2). Each band starts with an / = level, 
(b) The single-particle levels for Cgo after the icosahedral splittings are added to the spectra of (a). 
The tenfold degenerate HOMO (h„) and the sixfold degenerate LUMOl (tin) and LUM02 (tig) 
are denoted; they originate from the spherical / = 5 and / = 6 (tig) tt levels displayed in panel 
(a). For the a electrons, the icosahedral perturbation strongly splits the I = 9 level of panel (a). 
There result five sublevels which straddle the cr-electron gap as follows: two of them (the eightfold 
degenerate g„, and the tenfold degenerate h„) move down and are fully occupied resulting in a shell 
closure (180 a electrons in total). The remaining unoccupied levels, originating from the I = 9 a 
level, are sharply shifted upwards and acquire positive values. 



where * denotes complex conjugation. 

We take the radial parameters K;(r) as constants, and determine their value by adjusting 
the icosahedral single-particle spectra e-™ to reproduce the pseudopotential calculation of 
Ref. Q, which are in good agreement with experimental data. Our spectra without and 



with icosahedral splitting are shown in Fig. 2a and Fig. 2b, 



respective 
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y. We find that a close 



261] is achieved when the 



reproduction of the results of ah initio LDA calculations 
Wigner-Seitz radius for the jellium background is 1.23 a.u. The shell corrections including 
the icosahedral splittings are calculated using the icosahedral single-particle energies ef" in 
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eq. ([5]). The average quantities (p and V) are maintained as those specified through the ETF 
variation with the spherically symmetric profile of eq. (jlj). This is because the first-order 
correction to the total energy (resulting from the icosahedral perturbation) vanishes, since 
the integral over the sphere of a spherical harmonic Y"/" (/ > 0) vanishes. 

RESULTS 

Ionization potentials and electron affinities 

Having specified the appropriate Wigner-Seitz radius and the parameters ni of the 
icosahedral crystal field through a comparison with the pseudopotential LDA calculations 
for the neutral Ceo, we calculate the total energies of the cationic and anionic species by 
allowing for a change in the total electronic charge, namely by imposing the constraint 

4vr J pirydr = 240 ± a; , (9) 

where p(r) is given by eq. (jlj). The shell-corrected and icosahedrally perturbed first and 
higher ionization potentials are defined as the difference of the ground-state shell- 
corrected total energies as follows: 

4™ = E'J^°{Ne = 2m-x;Z = 240) - ^r(^e = 240 - x + 1; Z = 240) , (10) 

where Ne is the number of electrons in the system and x is the number of excess charges on 
the fuUerenes (for the excess charge, we will find convenient to use two different notations 
X and z related clS X — I Z \. A negative value of z corresponds to positive excess charges). 
Z = 240 denotes the total positive charge of the jellium background. 

The shell-corrected and icosahedrally perturbed first and higher electron affinities A'^° 
are similarly defined as 

= E'^°{Ne = 240 + X - 1; Z = 240) - ^^^(^e = 240 + x; Z = 240) . (11) 

We have also calculated the corresponding average quantities and A^, which result from 
the ETF variation with spherical symmetry (that is without shell and icosahedral symmetry 
corrections). Their definition is the same as in Eqs. fjTOj) and fjTTj) . but with the index sh 
replaced by a tilde and the removal of the index ico. 
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TABLE I: ETF (spherically averaged, denoted by a tilde) and shell-corrected (denoted by a 
superscript ico to indicate that the icosahedral splittings of energy levels have been included). IPs 
and EAs of fullerenes Cg^. Energies in eV. = 1.23 a.u. 

7 rico 7 4 ico 



1 5.00 7.40 2.05 2.75 

2 7.98 10.31 -0.86 -0.09 

3 10.99 13.28 -3.75 -2.92 

4 14.03 16.25 -6.60 -5.70 

5 17.09 19.22 -9.41 -8.41 

6 20.18 22.20 -12.19 -11.06 

7 23.29 25.24 -14.94 -14.85 

8 26.42 28.31 -17.64 -17.24 

9 29.57 31.30 -20.31 -19.49 

10 32.73 34.39 -22.94 -21.39 

11 35.92 39.36 -25.53 -22.93 

12 39.12 42.51 -28.07 -23.85 

In our calculations of the charged fuUerene molecule, the value and the icosahedral 
splitting parameters {ki, see Eq. ([7]), and discussion below it) were taken as those which 
were determined by our calculations of the neutral molecule, discussed in the previous sec- 
tion. The parameters which specify the ETF electronic density (Eq. (jlj)) are optimized for 
the charged molecule, thus allowing for relaxation effects due to the excess charge. This 
procedure is motivated by results of previous electronic structure calculations for CgQ and 



^60 
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26|, which showed that the icosahedral spectrum of the neutral Ceo shifts almost 
rigidly upon charging of the molecule. 

Shell-corrected and ETF calculated values of ionization potentials (IPs) and electron 
affinities (EAs), for values of the excess charge up to 12 units, are summarized in Table 
I (for Ts = 1.23 a.u.). Our shell-corrected results for the first electron affinity anc 



ionization potential are in good agreement with the experimental values (2.75 eV 4l| and 



first 



7.54 eV 421], respectively). For the second affinity, we find a small negative value of —0.09 
eV indicating that Cgg is a metastable species (for lifetime estimates, see below). For the 
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second ionization potential our value of 10.31 eV is close to that determined experimentally 
by McElvany et al. (9.7 eV) ll|] using charge-transfer bracketing experiments, and most 
recently by Sai Baba et al. (10.3 eV) [l^ using electron impact ionization. There is a 
controversy about this value, since several other experiments have found values that vary 
from 8.5 eV to 12.25 eV (for a detailed review of related experimental measurements available 

n 

in the literature, cf. Ref. [Ij], and references therein). Our theoretical calculations support 
the measurements of McElvany et al. and Sai Baba et al.., epecially since another quantity - 
namely, the appearance energy of Cgd", that is, the sum 

of E'=i I'T - which was measured 
12| is in very good agreement with the value calculated by us. 



in the recent work of Ref. 



Indeed, the appearance energy of Cg^J" determined in Ref. 1^ by linear extrapolation of the 
ionization efficiency curves was found to be 33.2±1 eV . This value should be considered 
as the upper limit, due to a curvature at the foot of the ionization efficiency curve. Ref. 



12l | estimates that this curvature leads to an overestimation of 2 eV . Taking this correction 
into account, the resulting experimental value of 31.2±1 eV is in close agreement with our 
theoretical value of 30.99 eV . 

The IPs and EAs calculated by us are plotted against the excess charge x in Fig. 3. For the 
ETF results (see Fig. 3a) the dependence is linear to a remarkable degree. Small deviations 
from linearity, however, are discernible and are due to the variations of electronic spill-out 
with varying excess charge. An inspection of the shell-corrected results in Fig. 3b reveals 
that they also exhibit a linear relationship with x within each electronic shell. Discontinuities 
due to shell openings are clearly discernible [ 1 < x < 10 for IPs corresponds to removal of 
electrons from the same h„ shell, while 11 < a; < 12 corresponds to removal of an electron 
from the shell immediately below it (g^ shell); 1 < x < 6 for the EAs corresponds to adding 
an electron to the same ti„ shell, while 1 < x < 12 corresponds to an extra electron in the 
same ti^ shell. The gap between the ti„ and the ti^ shell is approximately 1.23 eV , while 
the gap between the h„ shell and the gg shell is 1.75 eV]. The bending of the shell-corrected 
affinities around x = 9 — 12 is due to the fact that the electron emission for such highly 
charged states of the free molecule is barrierless and thus the uncertainty of the LUMO 
energy is high (see below). In the case of intercalated or endohedral compounds, the LUMO 
is drastically more bounded, due to the attraction of the ionic intercalant, and one should 
expect a strong linear relation even for this range of excess charge. 
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FIG. 3: (a) Extended Thomas-Fermi LDA Ionization potentials [1^ and electron affinities {—A^ 
as a function of the excess charge x. 

(b) Shell-corrected icosahedral IPs (I^™) and —EAs (— ^^™) as a function of the excess charge x. 
Charging energies and capacitance 

Fig. 4a shows that the variation of the total ETF-LDA energy difference (appearance 
energies) lS.E{z) = E{z) — E{0), as a function of excess charge z {\z\ = x), exhibits a 
parabolic behavior. The inset in Fig. 4a exhibiting the quantity 

E{z)-E{0) 



(12) 



plotted versus z (open squares), shows a straight line which crosses the zero energy line at 
z = 1. As a result the total ETF-LDA energy has the form, 



E{z) = E{0) + '^'J^'' - A,z 



(13) 



Equation (fT3|) indicates that fuUerenes behave on the average like a capacitor having a 
capacitance C (the second term on the rhs of eq. ([T3|) corresponds to the charging energy of 
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FIG. 4: (a) ETF-LDA total energy differences (appearance energies) AE{z) = E{z) — E{0) as a 
function of the excess charge z (z < corresponds to positive excess charge). 
Inset: The ETF function g{z) (open squares), and the shell-corrected function g]^(z) (filled circles). 
For z > 1 the two functions are almost identical. 

(b) magnification of the appearance-energy curves for the region —2 < z < 4. Filled circles: 
shell-corrected icosahedral values [AEl^{z) = El^{z) — El^°{0)]. Open squares: ETF-LDA values 
[AE{z) = E{z) - mi 

a classical capacitor, corrected for the self- interaction of the excess charge 0, [sl). We remark 
that regarding the system as a classical conductor, where the excess charge accumulates on 
the outer surface, yields a value of C = 8.32 a.u. (that is the outer radius of the jellium 
shell). Naturally, the ETF calculated value for C is somewhat larger because of the quantal 
spill-out of the electronic charge density. Indeed, from the slope of g{z) we determine 46] 
C = 8.84 a.u. 

A similar plot of the shell-corrected and icosahedrally modified energy differences 
AE'l'i^iz) = El'^iz) - ^^h°(0) is shown in Fig. 4b (in the range -2 < z < 4, filled cir- 
cles). The function gl^{z), defined as in eq. (fT2l) but with the shell-corrected quantities 
{AE^^{z) and v4'™), is included in the inset to Fig. 4a (filled circles). The shift discernible 
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between g^^{—l) and 5'sh°(l) is approximately 1.7 eV, and originates from the difference 
of shell effects on the IPs and EAs (see Table I). The effect of shell-closures for z = —10 
and z = 6 (which was discussed earlier in the context of Fig. 3b for higher IPs and EAs) 
is not discernible in the case of g^^{z) due to the scale of the inset. The segments of the 
curve g^^{z) in the inset of Fig. 4a, corresponding to positively {z < 0) and negatively 
{z > 0) charged states, are again well approximated by straight lines, whose slope is close 
to that found for g{z). Consequently, we may approximate the charging energy, including 
shell-effects, as follows, 

E'^i^) = ^r(o) + ^^^^^ - ^r^ , (14) 

for negatively charged states, and 

^r(^) = ^r(o) + ^^^^^ + I'r^ , (15) 

for positively charged states. Note that without shell-corrections (i.e., ETF) Ii — Ai = 
e^/C = 27.2/8.84 eV 3.1 eV, because of the symmetry of eq. ( fT3l) with respect to z, 
while the shell-corrected quantities are related as — A^-^° ^ e^/C + Ash, where the shell 
correction is A,,, ^ 1.55 eV (from TABLE I, Jj™ - ^ 4.65 eV). 

Expression ([T^ for the negatively charged states can be rearranged as follows (energies 
in units of eV), 

El'^ix) - E'^°{0) = -2.99 + 1.54(a; - 1.39)^ , (16) 



in close agreement with the all-electron LDA result of Ref . [27 1 . 

Equations ( JT^ and ( JT5l) can be used to provide simple analytical approximations for the 
higher IPs and EAs (see the definition in eqs. ( ITOi) and ( ITT]) ). Explicitly written, A'^° = 
E'^{x - 1) - E'^{x) = Af" -{x- l)eVC and = /f° + (x - l)e^/C. Such expressions 
have been used previously [29] with an assumed value for C ~ 6.7 a.u. (i.e., the radius of 
the Cqo molecule, as determined by the distance of carbon nuclei from the center of the 
molecule), which is appreciably smaller than the value obtained by us (C = 8.84 a.u., see 
above) via a microscopic calculation. Consequently, using the above expression with our 
calculated value for Af° = 2.75 eV (see TABLE I), we obtain an approximate value of 
A2° = —0.35 eV (compared to the microscopically calculated value of —0.09 eV given in 



TABLE I, and —0.11 eV obtained by Ref. |27|) — indicating metastability of Cqq — while 



14 



employing an experimental value for Af° = 2.74 eV, a value of Ai^" = 0.68 eV was calculated 



in Ref. 



29|. 



Concerning the cations, our expression fllSI) with a calculated /{™ = 7.40 eV (see TABLE 
1) and C = 8.84 a.u. yields approximate values 18.5 eV and 31.5 eV for the appearance 
energies of CgQ and CgQ (compared to the microscopic calculated values of 17.71 eV and 
30.99 eV, respectively, extracted from TABLE I, and 18.6 eV for the former obtained in Ref. 



251]). Employing an experimental value for = 7.54 eV, corresponding values of 19.20 



eV and 34.96 eV were calculated in Ref. 



. As discussed in Ref. [l^, these last values 



are rather high, and the origin of the discrepancy may be traced to the small value of the 
capacitance which was used in obtai„,„g these estimates Ref. 0. 

A negative value of the second affinity indicates that CgQ is unstable against electron 



autodetachment. In this context, we note that the doubly negatively charged molecule C 



has been observed in the gas phase and is believed to be a long-lived metastable species [ig, 



2- 



17l |. Indeed, as we discuss in the next section, the small LDA values of found by us and 
by Ref. yield lifetimes which are much longer than those estimated by a pseudopotential- 
like Hartree-Fock model calculation [16| , where a value of ~ 1 yus was estimated. 



Lifetimes of metastable anions, CgQ 

The second and higher electron affinities of Cgo were found to be negative, which implies 
that the anions Cqq with x > 2 are not stable species, and can lower their energy by emitting 
an electron. However, unless the number of excess electrons is large enough, the emission 
of an excess electron involves tunneling through a barrier. Consequently, the moderately 
charged anionic fullerenes can be described as metastable species possessing a decay lifetime. 

To calculate the lifetime for electron autodetachmant, it is necessary to determine the 
proper potential that the emitted electron sees as it leaves the molecule. The process is 
analogous to alpha-particle radioactivity of atomic nuclei. The emitted electron will have a 
final kinetic energy equal to the negative of the corresponding higher EA. We estimate the 
lifetime of the decay process by using the WKB method, in the spirit of the theory of alpha- 
particle radioactivity, which has established that the main factor in estimating lifetimes is 
the relation of the kinetic energy of the emitted particle to the Coulombic tail, and not the 
details of the many-body problem in the immediate vicinity of the parent nucleus. 
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Essential in this approach is the determination of an appropriate single-particle potential 
that describes the transmission barrier. It is well known that the LDA potential posseses 
the wrong tail, since it allows for the electron to spuriously interact with itself. A more ap- 
propriate potential would be one produced by the Self-Interaction Correction (SIC) method 



of Ref. 



431]. This potential has the correct Coulombic tail, but in the case of the fullerenes 



presents another drawback, na mely Koopman' s theorem is not satisfied to an extent ad- 



equate for calculating lifetimes. 



471 ] In this context, we note that Koopman' s theorem is 



known to be poorly satisfied for the case of fullerenes even in Hartree-Fock calculations 30]. 
Therefore, the HOMO corresponding to the emitted electron, calculated as described above, 
cannot be used in the WKB tunneling calculation. 

Since the final energy of the ejected electron equals the negative of the value of the 
electron affinity, we seek a potential that, together with the icosahedral perturbation, yields 
a HOMO level in Cg^ with energy —A'^°. We construct this potential through a self- 
interaction correction to the LDA potential as follows, 

^VKB = ^lda[p] - Vh[^] - Kc[e|-] , (17) 

where the parameter ^ is adjusted so that the HOMO level of Cqq equals — ^4^^^°. In the above 
expression, the second term on the rhs is an average self-interaction Hartree correction which 
ensures a proper long-range behavior of the potential (i.e., correct Coulomb tail), and the 
third term is a correction to the short-range exchange- correlation. 

For the cases of CgQ and Cgg such potentials are plotted in Fig. 5. We observe that they 
have the correct Coulombic tail, namely a tail corresponding to one electron for CgQ and to 
two electrons for CgQ . The actual barrier, however, through which the electron tunnels is 
the sum of the Coulombic barrier plus the contribution of the centrifugal barrier. As seen 
from Fig. 5, the latter is significant, since the HOMO in the fullerenes possesses a rather 
high angular momentum (/ = 5), while being confined in a small volume. 



44| . we estimate for CgQ a macroscopic half-life of 



Using the WKB approximation 
4 X 10^ years, while for CgQ we estimate a very short half-life of 2.4 x 10^^^ s. Both these 
estimates are in correspondence with observations. Indeed, CgQ has not been observed as a 
free molecule, while the free CgQ has been observed to be long lived 16, l3] and was detected 
even 5 min after its production through laser vaporization 17]. 

We note that the WKB lifetimes calculated for tunneling through Coulombic barriers are 
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FIG. 5: WKB effective barriers used to estimate lifetimes for Cgg (a) and CgQ (b). Dashed lines 
correspond to barriers due solely to Coulombic repulsion and solid lines to total barriers after 
adding the centrifugal components. The thick horizontal solid lines correspond to the negative 
of the associated electron affinities (a) and (b). In the case of CgQ [panel (a)], the 
horizontal solid line at —^2™ = 0.09 eV crosses the total barrier at an inside point Ri = 9.3 a.u. 
and again at a distance very far from the center of the fullerene molecule, namely at an outer 
point R2 = -eV^r = 27.2/0.09 a.u. = 302.2 a.u. This large value of R2, combined with the large 
centrifugal barrier, yields a macroscopic lifetime for the metastable CgQ (see text for details). 

very sensitive to the final energy of the emitted particle and can vary by many orders of 
magnitude as a result of small changes in this energy, a feature well known from the alpha 



radioctivity of nuclei [4^. Since the second electron affinity of Cgo is small, effects due to 
geometrical relaxation and spin polarization can influence its value and, consequently, the 
estimated lifetime. Nevertheless, as shown in Ref. 271, inclusion of such corrections yields 

n 

again a negative second affinity, [48( but of somewhat smaller magnitude, resulting in an 



28!1, the stabilization effect of the 



even longer lifetime. Furthermore, as discussed in Ref. 
Jahn- Teller relaxation for the singly-charged ion is only of the order of 0.03 - 0.05 eV. Since 
this effect is expected to be largest for singly- charged species, CgQ is not expected to be 



17 



influenced by it [27 1. 

On the other hand, generahzed exchange-correlation functional with gradient corrections 
yield slightly larger values for the second electron affinity. For example, using exchange- 
correlation gradient corrections, Ref. found = —0.3 eV, which is higher (in absolute 
magnitude) than the value obtained without such corrections. This value of —0.3 eV leads 
to a much smaller lifetime than the several million of years that correspond to the value of 
—0.09 eV calculated by us. Indeed, using the barrier displayed in Fig. 5a, we estimate a 



lifetime for CgQ of approx. 0.37 s, when A. 



ICO 

2 



-0.3 eV. We stress, however, that even this 



lower-limit value still corresponds to macroscopic times and is 5 orders of magnitude larger 



than the estimate of Ref. 



16| . which found a lifetime of 1 fis for A 



ICO 

2 



-0.3 eV, since it 



omitted the large centrifugal barrier. Indeed, when we omit the centrifugal barrier, we find 
a lifetime estimate of 1.4 fis, when = —0.3 eV. 

SUMMARY 



We developed a new LDA method for investigation of complex carbon clusters, and, as an 
illustration, applied it to multiply charged fuUerenes Cg^. The main elements of this method 
are: (i) Use of a stabilized jellium (structureless pseudopotential) approximation [20] for the 
ionic background, instead of the standard jellium approximation; (ii) Use of a recently 
introduced, shell-correction method 0, (iii) Inclusion of the effect of the icosahedral 
symmetry via a perturbative treatment. 

For the neutral Ceo molecule, the results obtained by our method are in good agreement 



with experimental o 



calculations 



24 



25 



Dservations, as well as with previous self-consistent electronic structure 



26| 



For the multiply charged fullerenes, charging energies were microscopically calculated for 
up to a; = 12 excess charges. A semiclassical interpretation of these results was developed, 
which viewed the fullerenes as Coulomb islands [2, ISj possessing a classical capacitance. A 
value of 8.84 a.u. was found for this capacitance, in contrast to the much smaller value of 6.7 
a.u. inferred previously from the particle-on-a-sphere model [29|. The calculated values for 
the first ionization potential and the first electron affinity agree well with the experimental 
ones. For the second and third ionization potentials, there exist substantial discrepancies in 
the experimental measurements 12| . Our calculations, reflecting the more accurate value of 



18 



the capacitance of Ceo, support the results from charge transfer bracketing experiments 11 1 



12|. 



and from direct ionization experiments by electron impact 

Employing the analogies with the case of alpha-particle radioactivity in nuclei, we found 
that the doubly c harg ed negative ion is a very long-lived metastable species, in agreement 
with observations 



16 



17l |. It decays through spontaneous electron emission, and its macro- 



scopic lifetime is the result of a superposition of a Coulombic barrier and a large centrifugal 
barrier, through which the emitted electron tunnels. 

Acknowlegdments 



This research is supported by the U.S. Department of Energy, Grant No. AG05- 
^6ER45234. 



[1] C. Coulomb, Memoires de V Academic, (1786) p. 67ff; ibid. (1787) p. 452; M. Faraday, 

Experimental research in electricity, (1839) (Dover, New York, 1965), Vok I, p. 360ff. 
[2] M.A. Kastner, Artificial atoms, in Physics Today, voL 46 (1993) 24. 
[3] Single Charge Tunneling, eds. H. Grabert and M.H. Devoret, (Plenum, New York, 1992). 
[4] C. Yannouleas and U. Landman, Chem. Phys. Lett. 210 (1993) 437. 
[5] C. Yannouleas and U. Landman, Phys. Rev. B 48 (1993) 8376. 

[6] For experimental studies of charged cationic clusters, see C. Brechignac, Ph. Cahuzac, F. 

earlier, and M. de Frutos, Phys. Rev. Lett. 64 (1990) 2893; C. Brechignac, Ph. Cahuzac, F. 

earlier, J. Leygnier, and A. Sarfati, Phys. Rev. B 44 (1991) 11 386; T.P. Martin et al. Chem. 

Phys. Lett. 196 (1992) 113, and references therein. 
[7] For a recent review of singly charged anionic metal clusters, see J. G. Eaton, L. H. Kidder, 

H. W. Sarkas, K. M. McHugh, and K. H. Bowen, in Nuclear Physics Concepts in the Study 

of Atomic Cluster Physics, R. Schmidt et al Eds., (Springer, Berlin, 1992), Lecture Notes in 

Physics Vol. 404, p. 291. 
[8] J. V. Coe, G. H. Lee, J. G. Eaton, S. T. Arnold, H. W. Sarkas, K. H. Bowen, C. Ludewigt, 

H. Haberland, and D. R. Worsnop, J. Chem. Phys. 92 (1990) 3980. 
[9] R. N. Barnett, U. Landman, C. L. Cleveland, and J. Jortner, J. Chem. Phys. 88 (1988) 4421; 



19 



ibid. 88 (1988) 4429; H.-P. Kaukonen, R. N. Barnett, and U. Landman, J. Chem. Phys. 97 
(1992) 1365. 

[10] D. Scharf, J. Jortner, and U. Landman, J. Chem. Phys. 87 (1987) 2716. 

[11] S.W. McElvany M.M. Ross, and J.H. Callahan, Mater. Res. Soc. Symp. Proc. 206 (1991) 

697; Acc. Chem. Res., 25 (1992) 162. 
[12] M. Sai Baba, T.S. Lakshmi Narasimhan, R. Balasubramanian, and C.K. Mathews, Int. J. 

Mass Spectrom. Ion Processes 125 (1993) Rl. 
[13] R.K. Yoo, B. Ruscic, and J. Berkowitz, J. Chem. Phys. 96 (1992) 911. 

[14] C.W. Walter, Y.K. Bae, D.C. Lorcnts, and J.R. Peterson, Chem. Phys. Lett. 195 (1992) 543. 
[15] H. Steger, J. de Vries, B. Kamke, W. Kamke, and T. Drewello, Chem. Phys. Lett. 194 (1992) 
452. 

[16] R.L. Hettich, R.N. Compton, and R.H. Ritchie, Phys. Rev. Lett. 67 (1991) 1242. 

[17] P.A. Limbach, L. Schweikhard, K.A. Cowen, M.T. McDermott, A.G. Marshall, and J.V. Coe, 

J. Am. Chem. Soc. 113 (1991) 6795. 
[18] S.N. Schauer, P. Wilhams, and R.N. Compton, Phys. Rev. Lett. 65 (1990) 625. 
[19] W. Kohn and L. J. Sham, Phys. Rev. 140 (1965) A1133. 
[20] J.R Perdew, H.Q. Tran, and E.D. Smith, Phys. Rev. B 42 (1990) 11 627. 
[21] H. Kuhn, J. Chem. Phys. 17 (1949) 1198. 
[22] G.A. Gallup, Chem. Phys. Lett. 187 (1991) 187. 

[23] M.R. Savina, L.L. Lohr, and A.H. Francis, Chem. Phys. Lett. 205 (1993) 200. 

[24] N. Troullier and J.L. Martins, Phys. Rev. B 46 (1992) 1754; J.L. Martins, N. Troullier, and 

J.H. Weaver, Chem. Phys. Lett. 180 (1991) 457. 
[25] Arne Rosen and Bo Wastberg, J. Chem. Phys. 90 (1989) 2525; Bo Wastberg and Arne Rosen, 

Physica Scripta 44 (1991) 276. 
[26] L. Ye and A.J. Freeman, Chem. Phys. 160 (1992) 415. 
[27] M.R. Pederson and A.A. Quong, Phys. Rev. B 46 (1992) 13 584. 
[28] V. de Coulon, J.L. Martins, and F. Reuse, Phys. Rev. B 45 (1992) 13 671. 
[29] Y. Wang, D. Tomanek, G.F. Bertsch, and R.S. Ruoff, Phys. Rev. B 47 (1993) 6711. 
[30] J. Cioslowski and K. Raghavachari, J. Chem. Phys. 98 (1993) 8734. 
[31] J. Hutter and H.P. Liithi, Int. J. Quantum Chem. 46 (1993) 81. 
[32] J. Kohanoff, W. Andreoni, and M. Parrinello, Chem. Phys. Lett. 198 (1992) 472. 

20 



[33] K. Yabana and G.F. Bertsch, Jellium model for Cqo, preprint (1992). 

[34] N. Van Giai and E. Lipparini, Z. Phys. D 27 (1993) 193. 

[35] M.J. Puska and R.M. Nieminen, Phys. Rev. A 47 (1993) 1181. 

[36] C. Yannouleas, E. Vigezzi, and R.A. Broglia, Phys. Rev. B47 (1993) 9849; C. Yannouleas, 
Chem. Phys. Lett. 193 (1992) 587; C. Yannouleas, P. Jena, and S. N. Khanna, Phys. Rev. B 46 
(1992) 9751; C. Yannouleas and R. A. Brogha, Phys. Rev. A 44 (1991) 5793; C. Yannouleas, 
R.A. Broglia, and M. Brack, Phys. Rev. Lett. 63 (1989) 255. 

[37] B. Grammaticos, Z. Phys. A 305 (1982) 257. 

[38] J. Harris, Phys. Rev. B 31 (1985) 1770. 

[39] R.C. Haddon, L.E. Brus, and K. Raghavachari, Chem. Phys. Lett. 125 (1986) 459. 
[40] M. Gerloch and R.C. Slade, Ligand field parameters (Cambridge Univ. Press, London, 1973). 
[41] S.H. Yang, C.L. Pettiette, J. Conceico, O. Chesnovsky, and R.E. Smalley, Chem. Phys. Lett. 
139 (1987) 233. 

[42] I.V. Hertel, H. Steger, J. de Vries, B. Weisser, C. Menzel, B. Kamke, and W. Kamke, Phys. 

Rev. Lett. 68 (1992) 784. 
[43] J. P. Perdew and A. Zunger, Phys. Rev. B 23 (1981) 5048. 

[44] A. I. Baz', Y. B. Zel'dovich, and A. M. Perelomov, Scattering, reactions, and decay in non- 
relativistic quantum mechanics (Israel Program for Scientific Translations Ltd., Jerusalem, 
1969). 

[45] Other aspects difi^erentiating our study from previous ones are: Consideration of the icosahe- 



dral perturbation, unlike Refs. 



34j, 



351 ]: consideration of 240 active electrons (as is the case 



2M), unlike Ref. 



331 ] which 

. 3 



with ah initio LDA pseudopotential claculations, see i.e., Ref. 

considered 360 active electrons, Ref. [34] which considered 60 active electrons, and Ref. 
which considered 250 active electrons. 

[46] Due to the changing spill-out with excess charge z, the capacitance should be written as 
C + 5{z). For our purposes here the small correction 5{z) can be neglected. 

[47] For certain systems, such as for example sodium clusters, an orbitally-averaged-like SIC treat- 
ment yielded highest-occupied-molecular-orbital (HOMO) energies for anions in adequate 



agreement with the calculated e. 
[48] The sign conventions in Ref. 



ectron affinities (see, Refs. 



271 ] are the opposite of ours. 



21 



